UNIVERSITY  OF  SHEFFIELD 


CO 

00 

00 


CM 


< 

I 

Q 

< 


department  of 
chemical 
engineering 
and 

fuel  technology 


FUNDAMENTAL  STUDY  OF  THREE  DIMENSIONAL 
TWO  PHASE  FLOW  IN  COMBUSTION 
SYSTEMS 

J.  SWITHENBANK,  B.C.R.  EWAN, 

F.  BOYSAN,  W.H.  AYERS 
S.V.  MALEBRAN 


FINAL  REPORT 


JUNE  1983 


* 


t.  REPORT  NUMBER' 


2.  GOVT  ACCESSION  NO. 


3.  RECIPIENT'S  CATALOG  NUMBER 


4.  TITLE  (and  Submit) 

FUNDAMENTAL  STUDY  OF  THREE  DIMENSIONAL 
TWO  PHASE  FLOW  IN  COMBUSTION  SYSTEMS 


S.  TYPE  OF  REPORT  4  PERIOD  COVERED 

FINAL  REPORT 

TO  JUNE  1983 _ 


6.  PERFORMING  O^G.  REPORT  NUMBER 


AUTHORf*; 


a.  CONTRACT  OR  GRANT  NUMBERfft) 


PROF.  J.  SWITHENBANK,  B.C.R.EWAN, 

F.BOYSAN,  W.H. AYERS 

AFOSR  -  82  -  0272 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Dept,  of  Chemical  Eng.  &  Fuel  Technology 
University  of  Sheffield,  Mappin  Street, 

Sheffield  SI  3JD,  England 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 

2308  -  A2 

It.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

12..  REPORT  DATE 

AFOSR/NA,  Building  410,  Bolling  Air  Force  Base 
D.C.  20332  .. 

20th  Dec.  1983 

13.  NUMBER  OF  PAGES 

78 

14.  MONITORING  AGENCY  NAME  4  ADDRESSf/l  different  from  Controlling  Office) 

EOARD,  223  Old  Marleybone  Road 

London  NW1  5TH,  England 

15.  SECURITY  CLASS,  (of  thle  report) 

15«.  oeclassification/oowngraoing 

SCHEDULE 

I  is.  DISTRIBUTION  STATEMENT  (of  thlB  Rmport)  | 

CLEARED  FOR  PUBLIC  RELEASE 


*7.  DISTRIBUTION  STATEMENT  (of  thm  abetted  entered  In  Block  30,  ff  different  from  Rmport) 


It.  SUPPLEMENTARY  NOTES 


19-  KEY  WORDS  (Continue  on  rararaa  alda  II  nacaaaary  and  Idantlty  by  block  number) 


COMBUSTION  MODELLING,  TURBULENCE  MODELLING ,  DUMP  COMBUSTOR 
SWIRLING  FLOW 


20.  ABSTRACT  (Continue  on  rmrormm  mldm  1/  nmc mmmmry  and  Identify  by  btock  numbmr) 

Details  are  presented  on  the  latest  developments  in  the  mathematical  modell¬ 
ing  of  turbulence  with  particular  reference  to  the  pressure  strain  transport 
term.  Comparisons  are  made  of  the  radial  variation  of  normal  and  shear  stress 
with  published  data  for  a  round  jet  and  for  decay  of  turbulence  for  selected 
flow  fields.  Work  is  also  reported  on  the  measurement  and  calculation  of  flow 
fields  inside  a  dump  combustor  using  swirl  and  baffle  stabilisation.  The  value 
of  the  different  turbulence  approximations  in  predicting  the  flow  are  discusse 


PD  I  JAN *73  1473  COITION  OF  I  NOV  «»  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Enlaced) 


Grant  Number: 


AFOSR  -  82  -  0272 


FUNDAMENTAL  STUDY  OF  THREE  DIMENSIONAL  TWO  PHASE  FLOW  IN 
COMBUSTION  SYSTEMS 


Prof.  J.  Swithenbank 

Department  of  Chemical  Engineering 

&  Fuel  Technology 

University  of  Sheffield 

England 


20th  December  1983 


Final  Report,  for  the  period  up  to  June  1983 


Approved  for  public  release 
distribution  unlimited 


Prepared  for 


AFOSR/NA 
Building  410 
Bolling  AFB 
D.C. 20332 


EUROPEAN  OFFICE  OF  AEROSPACE  RESEARCH  AND  DEVELOPMENT 
London,  England 


r 


Accession  For 

NTIS  GRAAI  W 

DTIC  TAB  fjj 

Unannounced  □ 

Justification _ 

By__ - 

Distribution/ _ 

Availability  Codes 
Avail  and/or 
Dist  Special 


o  * 


Introduction 


The  overall  objective  of  this  program  has  been  to  carry  out  fundamental 
research  to  support  a  continuing  design  program  for  combustors  of  interest 
to  USAF. 

* 

This  goal  is  one  which  has  many  contributing  elements,  and  of  necessity, 
the  emphasis  of  the  work  in  any  one  period  has  moved  to  address  the  most 
relevant  aspects  of  the  problem.  In  broad  terms  these  have  been  seen  to 
be  :** 

1.  The  refinement  of  the  individual  models  which  contribute  to  the 
overall  finite  difference  combustion  algorithm.  These  include  part¬ 
icularly  kinetics,  fuel  preparation  and  evaporation  and  the  radiation 
flux  contribution  to  enthalpy. 

2.  The  fundamental  analysis  of  an  adequate  turbulence  representation. 

3.  Research  into  efficient  iterative  procedures. 

4.  To  study  the  possibilities  for  interactive  software  suitable  for 
non-expert  usage. 

5.  Carrying  out  of  measurements  to  validate  models. 

The  following  review  describes  the  progress  made  in  recent  years  in  each 
of  these  areas  and  emphasises  the  most  important  areas  for  future  work. 

Review  of  Recent  Progress 

Combustion  Modelling 

Fuel  preparation,  evaporation  and  mixing  are  of  central  importance  to 
any  combustion  model. 

The  evaporation  problem  was  first  approached  from  the  viewpoint  of  indiv¬ 
idual  droplets  by  considering  the  trajectories  of  droplets  through  a  given 
combustor  flow  field  based  on  suitable  droplet  drag  and  evaporation  laws 
applied  to  each  finite  difference  cell. 

This  was  followed  by  an  extension  to  three  dimensions  and  the  integration 
of  both  Lagrangian  and  Eulerian  solution  procedures  for  droplets  and  flow 
fields  into  an  iterative  scheme.  This  enabled  one  to  observe  the  fate 
of  individual  droplet  size  groups  for  a  particular  combustor  geometry, 
and  results  were  presented  for  hot  and  cold  flow  fields  for  a  gas 
turbine  can,  demonstrating  different  droplet  trajectory  patterns  for 
different  spray  cone  angles. 

It  thus  became  possible  to  match  nozzle  type  with  combustors  since  un¬ 
evaporated  droplets  escaping  from  the  combustor  could  be  monitored  for 
unsuitable  nozzles.  This  has  been  seen  as  an  important  achievement  and 


has  generated  interest  from  several  quarters  including  Parker-Hannefin, 
Lycoming  and  Pratt  &  Whitney. 

Following  the  earlier  developments,  the  interaction  of  droplets  and  flow 
is  now  fully  coupled  in  both  directions.  In  addition,  the  influence  of 
turbulence  on  droplet  motion  is  incorporated  by  applying  a  randomised 
turbulent  contribution  to  the  mean  flow  velocity  in  each  cell,  based  on 
local  turbulence  intensity  and  a  random  variable.  This  is  important  in 
representing  the  spread  of  droplets  in  low  mixing  regions. 

An  optimised  iterative  scheme  for  the  combustion  case  is  summarised  as 
follows  - 

i)  The  u,v,w  momentum  equations  are  each  solved  in  turn  using  guessed 
pressures 

ii)  Since  the  velocities  do  not  satisfy  the  mass  continuity  equation 
locally,  a  "Poisson  type”  equation  is  derived  from  the  continuity 
equation  and  the  three  linearised  momentum  equations.  This  pressure 
correction  equation  is  then  solved  to  obtain  the  necessary  correct¬ 
ions  to  the  pressure  field  and  corresponding  adjustments  to  the  vel¬ 
ocity  components  are  made. 

iii)  The  k  and  £  equations  are  solved  using  the  updated  velocity  field 
to  obtain  the  distribution  of  effective  viscosity. 

iv)  Any  auxiliary  equations  (e.g.  enthalpy,  species,  conservation,  rad¬ 
iation  or  turbulence  properties)  are  solved  using  the  previously 
updated  values  of  the  other  variables. 

v)  Where  interphase  coupling  is  to  be  included,  a  solution  of  the 
equations  of  motion  of  the  droplets  is  obtained  at  intervals  and 

used  to  augment  the  source  terms  of  the  appropriate  gas  flow  equations 
These  steps  can  be  continued  until  the  error  mentioned  in  (ii)  has  dec¬ 
reased  to  the  required  value. 

Turbulence  Modelling 

It  has  been  seen  as  important  that  the  approach  to  the  turbulence  repres¬ 
entation  be  of  a  fundamental  nature  in  order  that  the  model  be  applicable 
to  a  wide  range  of  geometries.  The  problem  revolves  around  the  represent¬ 
ation  of  the  Reynolds  stress  components  and  the  various  options  which  are 
commonly  available  in  turbulence  modelling  are  shown  in  Figure  1. 

These  broadly  divide  between  effective  viscosity  models  and  fundamental 
models.  In  the  former,  the  individual  stresses  are  related  to  an  isotropic 
turbulence  viscosity  field.  This  in  turn  is  related  to  individual  turbul¬ 
ence  parameters  such  as  kinetic  energy  k  and  dissipation  rate  £  in  the 


k  -  e  model,  which  are  solved  from  their  respective  transport  equations. 

For  low  swirl  cases,  the  k  -  e  model  has  been  modified  in  the  past  by 
reducing  the  effective  viscosity  for  stresses  involving  the  swirl  component. 
This  was  found  however  to  have  only  limited  success  and  the  more  fundamental 
option  detailed  in  Figure  1  has  been  pursued  at  length. 

Much  of  the  development  work  on  turbulence  has  been  made  with  reference  to 
simple  swirling  flows  such  as  cyclone  flows.  One  of  the  most  striking 
demonstrations  of  the  shortcomings  of  the  k  -  e  model  was  shown  in  these 
flows  in  which  the  reverse  axial  flow  on  the  axis  is  not  represented.  It 
was  also  shown  that  by  contrast,  the  algebraic  model  represented  quite 
closely  both  the  mean  flow  contour  and  the  turbulence  throughout  the  cyclone. 
This  geometry  is  a  fairly  simple  one  mathematically  since  the  axial  deriv¬ 
atives  are  small  and  in  some  cases  can  be  neglected.  For  more  complicated 
flows,  such  as  dump  combustors  with  swirl,  more  solution  space  and  time  is 
required.  The  above  results  however  ,  have  been  very  rewarding  for  the 
considerable  effort  expended  and  have  indicated  strongly  the  direction  in 
which  combustion  model  development  must  proceed  if  the  representation  of 
swirl  geometries  mathematically  is  to  be  adequate. 

Interactive  Software 

An  important  part  of  any  intricate  design  problem  is  the  need  to  be  able  to 
interact  with  the  software  both  for  problem  definition  and  during  calc¬ 
ulation. 

A  parallel  development  procedure  has  therefore  been  undertaken  to  ensure 
that  problems  can  be  set  up  in  a  matter  of  hours  and  that  there  is  easy 
access  during  runtime  to  monitor  and  change  variables. 

This  is  best  illustrated  with  refereiice  to  Figure  2. 

Combustor  geometries  are  not  defined  analytically  but  rather  by  the  choice 
of  definition  of  each  of  the  cells  in  the  finite  difference  grid.  This  is 
a  simple  matter  of  assigning  letters  to  each  cell  and  some  options  include 
W  for  wall,  I  for  inlet,  0  for  outlet  etc.  An  arbitrary  geometry  can  thus 
be  built  up.  The  program  also  provides  a  list  of  options  at  each  stage  as 
well  as  a  HELP  command,  which  can  be  used  at  any  point,  to  help  the  inexper¬ 
ienced  user. 

Experimental  Measurements 

These  have  mainly  been  concerned  with  the  validation  of  particular  aspects 
of  model  refinement  and  are  complementary  to  data  used  from  the  literature. 
They  consist  of  the  measurement  of  temperature  profiles  to  compare  with 
the  predictions  of  the  combustion  algorithm,  residence  time  functions  and 
velocity  profiles  to  compare  with  the  predictions  of  the  mean  flow  and 
turbulence  models. 


In  the  residence  time  study,  a  short  pulse  of  mercury  vapour  (1  msec) 
was  used  as  tracer  in  a  gas  turbine  can  and  monitored  at  the  exit.  By 
different  choice  of  input  location,  the  different  stirred  and  plug  flow 
sections  could  be  studied.  In  itself,  this  is  a  quick  and  useful  method 
of  assessing  reactor  volumes  and  blow  off  performance. 

Results  however,  were  compared  with  mixing  predictions  for  the  k  -  e 
model  using  the  randomised  bubble  tracking  procedure  developed  for  partic¬ 
les  and  referred  to  above.  Of  the  five  input  locations  studied,  good 
agreement  was  obtained  for  the  downstream  positions  with  regard  to  si  ape 
and  position  of  the  response  function.  For  the  upstream  input  positions, 
good  agreement  with  shape  was  obtained,  indicating  that  the  mixing  volume 
is  well  represented  in  the  model.  A  slight  discrepancy  in  plug  flow  times 
for  upstream  input  was  thought  to  be  due  to  some  computational  trapped 
volume  on  the  symmetry  axis  of  the  combustor. 

The  method  has  proved  a  useful  one  in  enabling  a  comparison  to  be  made 
with  predictions  at  the  fundamental  mixing  level. 

The  measurement  of  flow  velocity  profiles  using  LDA  has  been  carried  out 
in  conjunction  with  the  turbulence  model  development  and  has  been  necessary 
since  little  experimental  data  existed  for  the  simple  swirl  flows  which 
the  development  program  required. 

For  this  purpose,  two  types  of  tangential  entry  cyclones  have  been  measured 
and  the  present  good  agreement,  which  the  algebraic  turbulence  model  has 
provided  has  already  been  reported. 


Progress  during  1982-83 

The  objectives  of  the  past  year  have  sought  to  address  the  particular  prob¬ 
lems  of  swirling  flows  with  reference  to  their  application  in  dump  comb¬ 
ustors  and  may  be  itemised  as  follows  : 

1)  Further  development  of  the  Sheffield  algebraic  turbulence  model  with 
particular  reference  to  the  pressure  strain  term,  the  role  of  which 

is  to  transfer  energy  between  the  different  fluctuating  velocity  comp¬ 
onents  tending  to  equalise  the  normal  Reynolds  stresses  and  decrease 
the  absolute  value  of  the  shear  stresses. 

2)  To  investigate  the  effects  of  different  types  of  swirl  on  the  velocity 
field  in  a  dump  combustor.  Swirl  types  woc^d  range  from  free  vortex 
type  to  solid  body. 

Investigations  would  include  - 

a)  application  of  a  number  of  flow  model  approximations  to  the  diff¬ 
erent  types  of  isothermal  flow  to  confirm  the  regimes  of  agreement 
and  difference. 


b)  application  of  laser  Doppler  anemometry  to  determine  the  flow  fields 
in  the  dun?  combustor  for  the  different  swirl  types. 

The  following  is  a  review  of  the  work  carried  out  in  each  of  these  areas 
and  begins  with  a  description  of  the  turbulence  modelling. 
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REPORT  ON  THE  MODELLING  OF  THE 
PRESSURE  INTERACTION  TERM 


Sergio  Vasquez  Malebran 

Department  of  Chemical  Engineering  and  Fuel 
Technology,  Sheffield  University. 


Abstract 

This  report  investigates  the  effects  of  the  non- 
linearities  on  the  quasi-isotropic  model  for  pressure-strain  in 
the  Reynolds  stress  equation.  In  particular  Rotta's  linear  return 
to  isotropy  model  has  been  replaced  by  an  expression  which  also 
contains  the  second  invariant  of  the  anisotropy  tensor  and  a 
quadratic  term  for  the  Reynolds  stresses.  The  constrlbutlon  of 
the  mean  strain  to  the  pressure  interactions  is  modelled 
according  to  Launder,  Reece  and  Rodi  (1975)  and  a  complete 
expression  of  the  fourth  order  tensor  is  considered.  A  simplified 
version  of  this  model  is  also  studied  where  the  parameter  c2  of 
the  standard  LRR  model  is  made  a  function  of  the  invariants  of 
the  anisotropy  tensor.  Constants  resulting  from  the  approximate 
expressions  are  adjusted  in  reference  to  homogeneous  turbulent 
flows  and  satisfactory  agreement  was  obtained  for  a  single  set  of 
constants . 

1 .  Introduction 

The  main  problem  in  predicting  turbulent  flows  is  to 
determine  the  Reynolds  stresses,  <UjUj>  which  appear  in  the 
equations  of  conservation  of  momentum.  Attempts  have  been  in  the 
past  to  evaluate  these  stresses  by  relating  them  to  known  or 
calculable  quantities  such  as  the  mean  velocity  gradients.  A 
simple  expression  of  this  type  which  is  applicable  to  situations 
in  which  there  is  a  single  dominant  direction  of  flow  is 


-<U] U2>=  Vi 


where  - —  is  the  only  non-zero  mean  velocity  gradient,  u  is  the 
OX  2 


3fr. 
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fluctuating  velocity  and  v*.  is  the  so  called  eddy-viscocity.  For 

I 

■ore  coaplex  flows  this  simple  relation  can  be  generalized  in  the 
following  fashion 


■<uiV  s  vt(i7 


2 

d»i  3 


where  k=— <u.u.>  is  the  turbulent  kinetic  energy.  The  second  term 

2  *  J 

on  the  right  hand  side  has  been  added  to  make  eq.(2)  contract 
properly.  Unlike  the  molecular  viscocity  ,  which  is  a  property  of 


the  fluid,  the  eddy  viseoeity  is  a  property  of  the  flow  and  many 
hypotheses  have  so  far  been  proposed  for  its  evaluation. 
According  to  the  Prandtl-Kolmogorov  hypothesis,  the  eddy 
viscocity  can  be  expressed  as 


vt  =  C  k*l  (3) 

where  k  denotes  the  turbulent  kinetic  energy,  1  a  length  scale 
proportional  to  that  of  the  energy  containing  motions  and  is  a 
constant  of  proportionality.  For  relatively  simple  flows  where 
only  one  Reynolds  stress  component  is  of  importance,  equations 
(2)  and  (3)  give  fairly  good  results.  For  flows  where  several 
components  of  the  Reynolds  stresses  are  of  Importance, 
experiments  have  shown  that  can  be  strongly  direction 
sensitive.  Furthemore,  the  empirical  coefficient  can  no  longer 
be  regarded  as  a  constant.  Thus  ,  the  turbulence  models  based  on 


i;  J 


the  eddy  viscocity  concept  which  is  expressed  in  equation  (2)  and 
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(3)  do  not  produce  the  universality  required  for  the  calculation 

I 

of  complex  thr*ee  dimensional  flows. 

The  governing  equations  of  the  transport  of  the 
Reynolds  stresses  can  be  derived  from  the  Navier-Stokes  equations 
and  can  be  expressed  in  the  following  way: 


where  upper  and  lower  case  u*s  refer  to  mean  and  fluctuating 
velocity  components  respectively,  p  is  the  pressure,  x’s  are  the 
cartesian  space  coordinates,  v  and  p  are  the  kinematic  viscocity 
and  the  density  of  the  fluid. The  physical  meaning  of  these  terms 
is  similar  to  that  of  the  corresponding  terms  of  the  equation  for 
the  turbulent  kinetic  energy  obtained  by  contracting  (4);  the 
term  on  the  left  hand  side  represents  the  convection  of  the 
Reynolds  stresses  (C^j)  while  the  ones  on  the  right  hand  side  are 
the  production  term  (P^)  which  represents  the  exchange  of  energy 
from  the  mean  motion  to  the  turbulent  motion,  the  dissipation 
term  (Ajj)  which  expressothe  diffusive  transport  of  the  Reynolds 
stresses,  the  pressure  interaction  term  (4^)  which  represents 
the  correlation  between  velocity  fluctuation  and  the  pressure 
gradient,  and  finally  the  viscous  dissipation  term  (e^j)  *  The 
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last  three  termp  require  modelling  assumptions  in  order  to  relate 

) 

them  to  known  or  calculable  quantities.  The  modelling 
assumptions,  however,  are  much  less  restrictive  in  this  case  than 
they  are  in  the  case  of  the  general  eddy  viscocity  hypothesis  and 
the  models  based  on  the  above  transport  equation  are  likely  to 
possess  a  much  greater  potential  for  predicting  complex  flows 
successfully.  Equation  <*)  is  in  fact  a  set  of  differential 
equations  for  the  individual  Reynold  stress  components  and  the 
solution  of  this  set  of  partial  differential  equations,  which  of 
necessity  is  a  numerical  one,  can  be  very  expensive  in  terms  of 
computing  time. 

2«  Approximated  expressions 


2. 1  Dissipation  Term 

The  simplest  model  for  the  dissipation  term  in  eq.(4) 
is  given  by 


(5) 


where 


This  expression  holds  only  when  the  dissipative  motions  are 
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isotropic  and  should  be  applicable  to  those  flow  regions  where 

I 

the  local  turbulence  Reynolds  numbers  are  high.  Although 
experimental  studies  have  shown  that  tur* ulenee  does  not  remain 
locally  isotropic  in  the  presence  of  strong  strain  fields,  this 
expression  has  been  used  in  previous  modelling  efforts  with 
reasonable  success. 

2.2  The  Diffusion  Term 

Although  the  requirement  of  high  Reynolds  numbers 
permits  the  viscous  diffusion  term  in  eq.(H)  to  be  neglected,  it 
is  not  possible  to  dismiss  the  remaining  diffusion  terms  with 
confidence.  Hanjalic  and  Launder  P97?)  suggest  that  the 
•pressure  diffusion*  terms  may  also  be  neglected.  Although  their 
assumption  is  based  on  a  single  result,  in  the  absence  of  any 
other  firm  evidence,  this  seem* to  be  the  best  hypothesis,  thus 
only  the  triple  velocity  correlations  may  be  replaced  by  at: 
expression  containing  only  the  second  order  correlations  in  the 
following  way: 

, .  d<uiu.>  d<u.  u.  >  d<u.u,> 

-<u1u.u.>  =  c  -(<u  u  .>- — - -  +  <u,u..> — ; -  +  <u.u.> — r — (&) 

i  j  k  »c  i  j  dXj  J  i  axj  k  i  axj 

> 

where  c  is  a  constant.  A  simpler  gradient  diffusion  hypothesis 
s 

has  been  proposed  by  Daly  and  Harlow  (1970)  which  is  expressed 
as : 
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'  k  6<uiuJ> 

•<U.U.U.>  =  C  —  <u.u,> - *— 

1  i  k  s£  k  1  a* 


where  c  is  a  constant.  Launder  Reece  and  Rodl  (1975)  suggest 

9 

» 

that  the  constants  c  and  c  take  the  values  of  0.11  and  0.22 

9  9 

respectively. 


2.3  Pressure  Interaction  Term 

for  mod  ling  purposes  the  pressure  interaction  term  is 
often  split  into  a  sum  of  two  components 


Bu.  Bu 


-  (-(u.^-  +  u  ^£-)}  =  (£(•— ±  -  /!(_&_( pu)  +  -A-tpu  ))\  (8! 

'p  J8xa  iBXj  '  xp  BXj  bx±  f  'p  BXj  J  Bx j  i  • 


where  the  first  term  on  the  right  hand  side  is  the  pressure- 
strain  correlation  and  since  its  trace  vanishes  its  role  is  to 
transfer  energy  among  components  and  not  to  create  or  destroy  it. 
The  second  term  is  referred  as  the  'pressure  diffusion'  and  it  is 
usually  neglected  in  the  modelling  procedures. 

Following  Chou  (19^5),  the  explicit  appearance  of  the 
pressure  fluctuation  satisfies  the  following  Poisson's  equation: 


1  B2p 

pBXjBXj 

with  solution 


B2  Bu.BU 

- — r — (  u  u.  -  <u,  u.  >  )  ♦  2- — - — 

Bxia*k  1  k  1  k  axi8xk 


»<  »  r  't  l{ 


535 
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p(x) 


'  +  - 

1 


4« 


-ai-  (u<Uk 

5x,&x.  i  k 


<u1uk>) 


♦  2 


5u^i 

BxidXk! 


”!(X») 


dV 


X'  -  X 


♦  riJ  <10> 


where  the  integration  la  carried  out  over  the  whole  at*  -  space 
and  specified  at  the  point  (jr ' )  ,  and  is  a  surface  integral. 
The  pressure-strain  correlation  may  be  expressed  in  the  form: 


pSX, 


1 

4* 


u  ,du 

—  >  r-1  +  2  (  r“ 

®x,  bx  6x.  6x 

1  O  J 


4V°i 


du 


i 

ii. 


dV 


ia,B  °*jj  ji-  -  *i 


ru  (11> 


this  equation  shows  that  the  pressure-strain  originates  from  the 
Interaction  between  turbulence  components  (  (4*^),^)  and  from  the 
interaction  of  the  mean  rate  of  strain  with  turbulence  (  < 4*^  j )  •  2 > 
.  The  surface  integral  is  negligible  away  from  the  vecinity  of  a 
solid  boundary. 

Most  of  the  modelling  work  so  far  has  adopted  Rotta's 
(1951)  model  for  (  (♦1J)M) 


where  cj 

kinetic 

quotient 


)  =  -e.-{  <u.u.>  -  —6.  .)  (12) 

M  a  *  i  3  1 J 

is  a  constant,  k  and  e  are  the  time  averaged  turbulent 

energy  and  the  dissipation  energy  respectively,  the 
|( 

—  thus  represents  a  characteristic  decay  time  of 
e 


turbulence 
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The  aean  strain .  contribution  to  the  correlation  ■*■ -  is 

,  *xi 

very  often  modelled  starting  from  the  quasl-isotropic  model.  This 

was  first  published  Hanjalic  and  Launder  (1972)  and  by  Launder  et 

al  (1973).  Other  derivations  following  dlferents  paths  have  been 

proposed  by  Naot,  Shavlt  and  Volfshteln  (1972)  and  Lumley  (1973). 

Launder  Reece  and  Rodi  (1975)  extended  this  model  to  take  in 

account  'wall-echo'  effects.  Erdogan,  Boysan  and  Swithenbank 

(1980  )  assumed  a  constitutive  equation  for  this  term,  the 

pressure-strain  being  a  function  of  the  mean  rate  of  strain 

tensor  E  ,  the  mean  rotation  tensor  Q  and  the  Reynolds  stress 

tensor  a  and  a  similar  expression  to  that  obtained  by  the  quasl- 

isotropic  model  is  reached  with  three  Independent  parameters 

being  made  as  a  function  of  the  local  ratio  of  production  to 

.  P 

dissipation  —  .  Lin  and  Wolfhstein  (1979)  have  assumed  a  complete 
e 

different  approach,  arriving  at  a  more  complicated  form  but  is 
work  still  remains  to  be  tested  in  more  than  simple  flow 
situation 


3.  Non-linear  Model  for  the  Pressure-strain  term 


3«1  The  Quasy-Isotroplc  Non-linear  Model 

In  simulating  ^ ^ ^ ^  1 1  Rotta's  proposal  can  be  expanded 
incorporating  non-linear  terms  as  follows: 


!-£<c1«;n)«„lV  -f*u>-c’i^(<uiuk><"kuj>  -i<“rv<u,ur>6ij> 
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»  '.II 

where  fcj  ,  cj  and  e^  are  constanst  to  be  adjusted*  and  II 
denotes  the  turbulence  anisotropy 


11  *  S7<  <ulV  -  f6U><  <"lV  -  I5!!*  =  »lkbkj 


(13) 


In  simulating  (<>..)  following  LRR  (1975)  it  is  assumed 

1 J  i2 

that  this  correlation  may  be  aproxiaated  in  the  form: 


,  .  .  .  .1^1 

ij  *2  1Jbxb 


where 


„mi 

•u  * 


1  f  a8%“l  dV 

' 2*  I  -  ,j 


(14) 


and  the  5,  are  the  cartesians  components  of  the  position 
vector  jii*  -  xj  .  This  expression  is  exact  for  homogeneous  flows 
but  aproxlmately  true  in  more  general  flows. 

Rotta  (1951)  has  commented  that  the  fourth  order  tensor  aB* 
should  satisfy  the  following  kinematic  costraints: 


(i)  a; 


alB  -  aBi 

"lj  '  *J1 


(Symmetry 


(id  =o 


(Continuity 
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(Green’s  Theorem) 


The  most  general  expression  for  this  quasy-isotropic  model 
including  non-linear  terms  is: 


*1J  *  “<Vl>SlJ  ♦  P(<umul>6iJ  *  <U.UJ>611  *  <uiW)  *  <uiuj>6.l) 
tl)  [2]  [33  [4]  [53 

+  c2<ulUj>6ni  +  4  T,(6ml6iJ  4  6mj6n)) 

[63  [73  [83  [93 


,<UU.>  I  <UmU1>  <««*> 

*  a  — iT-<“l“j>  *  f  <— — ' <U1"J>  *  — iTJ-<0lUl>) 


[103 


[113 


[123 


*  ♦  V(<U»Uk><UkUl>6iJ  *  <“»“k><UKUJ>6ll 


k  k  i  mj  k 
[133  [143 


[153 


4  <uiukXujuk>6Bl  4  <uiukXu1uk>6|nj) 


[163 


[173 


»<ulV 

*  c2— T~<UjUk>6mi  + 
[183 


(15) 


<u  u  ><u  u  >  *  *  *  \  M, 

*  ”  (6»i6u  *  V)6n>)k 

k2  Bi  1J  [203  [213 


[193 


>_«  o,,  +  » 

1  1 J  [203 


5? 


where 
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a,  p.  a/t  p  ,  a  ,  p  ,  v,  r),  v  ,  n  ,  c2,  and  c2  are  constants, 


A  linear  expression  can  be  obtained  by  simply  putting 

tilt  •  n  k 

a  *  p  <  v  •  t)  ,  e2  ,  a  and  p  equals  to  zero  and  by  using 

the  constraints  one  obtains 


o<u  u->  +  5p<u  u  >  +  e2<u  u.>  +  k6  .(2p  +  v  ♦  4tj)  =  0  (ii) 


(Ba  +  i»p)<uinu1>  +  k  (  2c2  +  3v  +  2n)&ni  =  2<u0ui>  (ill) 


which  gives  4  equations  in  5  unknowns 


c  ♦  5p  ♦  c2  =0 
2B  ♦  v  ♦  4t)  =0 
3a  ♦  4p  =0 
2c2  +  3v  ♦  n  =0 


(I6.a) 


a  t  P  •  v  and  t)  can  be  determlneds  in  term  of  c2  as: 


a  =  y^(4e2  ♦  10)  p  =  — —(  3c2  +  2) 


n  =  ~M20c2  ♦  2) 
55 


v  =  -tt(50c2  +  4) 
55 


In  a  similar  way  also  the  coefficients  for  the  linear  +  non 
linear  model  can  be  obtained  as  follows. 


'•v**n*#*f  - 


J 
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■k  . 
*lk  =0  ! 


(Continuity) 


(o+  5P>  c2  ♦  2p’)<ujnu1>  ♦  k6Bl(2p  ♦  v  +  «n)  ♦  ... 


<UBUk>.  ..»»«*  »  »  (O'AXu.u^)  -  .. 

+  __<Uiuk>(a  +p  +a  +5p  +C2>  +  - LCS  P  *6.1(p"+v,*W ) 


giving 


a  *  5p  +  c2  +  2p  =  0 
2p  +  v  +  !»t)  a  0 


(16. b) 


•  »  «  n  t 

o  +p  ♦  a  +5p  ♦e2=0 


«  i  » 

P  ♦  v  +  Htj  =  0 


and  from  the  Green’s  Theorem  a®*  =  2<u  u  > 

kk  mi 


2<Umui>  =  <umui><3a  ♦  4p  +  2a’)  +  6  ,k(2c2  +  3v  ♦  2n)  ♦ 


<u  u,  > 
m  k 


t  «t  „  <u  u  ><u  u  > 

k  -<ukui>«2P  -  3a  ♦  *p  )  +  k6  .  (c2  ♦  3v’  ♦  2r, 


3a  +  <p  4  2a  =2 
2c2  4  3v  4  2t|  s  0 


(16. c) 
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I  If  M  . 

2p  *  3a  *  Mp  =  0 


tit 

c2  ♦  3v  +  2r\  -  0 


The  result  Is  a  system  of  8  equations  with  11  unknowns 


Refering  to  (15)  we  and  collecting  terms  to  get 


•  i 


( 4»)  » 2  =  ~^a  * 

p>pu  - 

ppkk6ij 

-  (  c  2 

+  p)Djj  +  (a  +P 

♦ 

(13 

133 

(23 

[63 

[53  [103 

tie] 

ft 

+  (a 

•  •  • 

ft 

♦  (P  ♦ 

c’>vij 

,  <Uiu  > 

-  *  „  * 
[11] 

(133 

r  is3 

[173 

tie] 

k  (v  ♦  q)  ♦  (v*  +  n* )<UpUq>^UqUp>  E,  ,  + 
*•*  +  K  (7)  [9)  1 1 93  (213  k2  13 

- 


p  Skk6ij 

( 1  *•  3 


numbers  in  square  brackets  underneath  the  terms  refer  to  those  in 
equation  (15);  with  and  Djj  haring  the  following 
meaning : 


dUj  dU 


■13  ' 


pl 3  -  -<<ulV^  *  <ujV^> 
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'  '  &Uk  aUk 

-<<uiuk>^  +  <u3uk>577> 


also  S..,  T..  and  V..  are  the  following  tensors 
iJ  iJ  ij 


1 

>  s  “(<1111.  Xu.  u  .>  ♦  <u  u  Xu  u  >)- - 

lj  k  lj  fflj  II  ax 

ID 


1  au,  au. 

T ,  .  =  “( <u.  u .  y  <u  u.>  — -  ♦  <u.u.Xu  u,> - ) 

ij  k  k  i  mkftx  k  j  m  k  gx 

m  in 


,  au,  au, 

VU  *  k(<ui°k><uiuk>i^  *  <ujuk><uiuk>^7) 


*nd  skk  *  Tkk  *  vkk 


It  can  be  noted  that  due  to  continuity  in  an 
Incompressible  fluid  the  constribntion  of  the  terms  [i|]t  [8], 
(163  and  [20]  is  zero;  the  term  [11]  is  incorporated  into  ((}>.,) 

1 J  *  i 

and  neglected  since  its  variation  is  similar  to  that  of  the 
second  Invariant.  Collecting  all  terms  we  obtain  the  final  form 
of  this  non-linear  model. 


1  ,<U  U  ><U  U  > 

4»i j  =  (<’ij>M  +  Al(PirXv6ui  ♦  <a2  ♦  A2  £-3IT-a-E-)kE 


i j-3  kk°i j 


U 


♦  A3(DiJ-3Dkk61^)  ♦  A|*(Si j“3Skk6iJ) 
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V*5CTirlT 


♦  A6^vj j~3Vkk61j) 


where 

i  i  i 

Aj  =  -(o+p)  A2  =  (v+q)  A2  =  (v  +n  ) 

*  i  *»  w 

A3  a  -(c2+p)  A4  S  (a  +p  )  A5  S  (a  +p  ) 


(16) 


vf  9 

A6  S  (p  +c2) 


By  comparing  equations  (I6.b)  and  (16. c)  with  equation  (16. a) 
from  the  linear  case  it  can  be  noted  that  only  the  constants  t' 
and  p'  are  absent  from  (16. a).  It  may  be  concluded,  although  not 
definettly,  that  a%  and  p*  are  the  principal  non-linear 
constributions  to  a^  and  so  we  can  neglect  a"  ,  p"  ,  q'  and  v’. 
Now  expressing  o,  p,  q,  and  v  in  terms  of  c2 ,  o'  and  P'  ie.: 


a  =  -j^-(*tc2+10(  1-a  '  )  +  8p  '  ) 

q  s  -^(20c2+(1-a’)+18p’) 
55 

therefore  the  coefficients  A^ 
become 

Aj  =  --^•(cj+Sd-o'j+Pp') 
*3  =  -~r( 8c2-2 ( 1-a ' ) -6p  * ) 


p  =  ~(3c2  +  2(1-a,)  +  6p') 

v  =  --^(50c2+4(1-o,)+12p') 

55 

A 2  .  A3  and  Aj,  of  equation  (16) 

A 2  =  --V(30c2-2(1-o,)-6p') 

55 

«  t 

A4  =  (c  +P  ) 


for  a  *  p  *  0  ,  Aj  ,  »  *3  and  A4  reduces  to  the  linear 

t 

{ 

quasi-isotropic  model  of  LRR. 

Equation  (16)  is  the  most  general  form  of  the  pressure- 
strain  correlation  and  an  optimization  for  all  the  constants  is  a 
difficult  task,  it  may  be  then  be  reasonable  to  look  for  a 
simplified  version.  In  the  next  section  a  new  approach  is 
considered  in  which  it  is  assumed  that  the  parameter  cj  of  the 
standard  LRR  Model  is  sensitive  to  the  anisotropy  of  the  flow. 

3.£  Pressure-strain  Model  with  £2.  £  function  of  the  invariants 

In  our  investigation  we  have  found  that  in  solving  the 
differential  equations  for  Harris  et  al.  flow  using  LRR  model,  a 
constant  value  for  C2  fail  to  predict  the  level  of  tne  shear 
stress  uv  and  consequently  giving  unsatisfactory  predictions  for 
the  other  Reynolds  stresses;  similarly  in  the  turbulence 
undergoing  axisymmetric  strain  the  linear  and  non-linear  quasi- 
laotropic  model  are  unable  to  reproduce  the  deflection  of  energy 
to  the  longitudinal  component.  It  is  possible  to  improve  this 
situation  without  affecting  the  predictions  for  the  other  flow 
considered,  by  assuming  that  C2  is  influenced  by  the  state  of 
anisotopy  of  the  flow  and  so  c2  is  taken  as  s  function  of  the 
second  and  third  invariant  of  the  anisotropy  tensor,  ie: 


c2  =  0.4  ♦  albikbkJ  +  a2b*vbnb 


ik  k j  ji 


where 
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<0iUJ>  2 

bjj  *  - —  ■*—  -  *2^^j  ,  ®i  and  82  are  constants. 


although  so  far  only  the  second  invariant  has  been  taken  into 
account . 


JLi  Homogeneous  Turbulent  Flows 

Homogeneous  turbulent  flows  are  physical  models  on 
which  the  turbulent  motion  is  simpler,  they  can  be  described  by 
coupled  ordinary  differential  equations.  Theses  flow 
configurations  are  the  "ideal"  cases  for  computational  modelling. 
Flow  of  this  class  may  be  devided  in  three  groups 


(1)  Homogeneous  Isotropic  turbulence.  This  is  the  simplest  type 
of  turbulence.  In  spite  of  the  fact  that  it  is  unrealistic  a 
knowledge  of  its  structure  and  behaviour  has  provided  the  basis 
for  a  better  understanding  of  non-isotropic  flows. 


(11)  A  second  physical  model  is  the  homogeneous  non-isotropic 
turbulence  with  no  mean  velocity  gradients.  In  this  situation 
there  is  a  strong  tendency  towards  isotropy  and  so  the  first  part 
of  the  pressure-strain  plays  a  significant  role.  For  this  and  the 
above  group  the  turbulent  kinetic  energy  is  always  decaying  since 
there  is  only  redistribution  and  dissipation  of  energy  (  we  hsve 
used  the  experimental  data  of  Uberoi  (1956)  and' Tuckers  and 


*5 


INTERNAL  REPORT 


-  18  - 


PRESSURE  HODELLING 


Reynolds  (1968)  ). 


(11)  A  third  physical  model  Is  the  homogeneous  non-ls  tropic 
turbulence  with  constant  mean  velocity  gradient.  In  this 
configuration  there  is  producction,  redistribution  and 
dissipation  of  energy.  We  may  subdivide  this  group  in  flows  of 


two  kind: 


Uniform  pure  strain  :  In  this  case  the  principal  axes  of  the 
Reynolds  stresses  are  the  same  as  those  of  the  mean  rate  of 
strain.  We  have  considered: 

Plane  strain  (  Townsend  (1956),  Tucker  and  Reynolds  (1968)  ) 
Axlsymetric  contraction  (  Uberoi  (1956)  Tan-atichat  (1980)  ) 


Uniform  rotational  strain  :  In  this  case  the  axes  of  the 
Reynolds  stresses  are  not  aligned  with  the  axes  of  the  mean  rate 
of  strain.  We  have  taken  for  this  kind  of  flow  the  experiments  of 
Champagne  et  al  (1970)  and  Harris  et  al .  (1977). 


Table  I  summarizes  the  flows  with  their  initial  conditions  which 
*  we  have  been  used  in  the  numerical  computation. 


5.  Numerical  Computation 
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Solving  the  Differential  Equations.  The  description  of 
these  turbulent  flows  was  done  by  solving  the  coupled 
differential  equations  for  the  Reynolds  stresses  with  their 


initial  conditions  and  the  dissipation  equation  : 


JL  -  uJL  ,  and  c  .  =  1.44  ,  c  ,  =  1.92 
dt  dx  ci  '  e2 

Programs  employing  Euler-First  order  method  were  used  except  for 
the  axlsymetric  contraction  where  a  variable-order,  variable-step 
Adams  method  (1976)  was  used.  The  presssure-strain  models  were 
Incorporated  in  these  programs  and  the  constants  were  adjusted  to 
obtain  satisfactory  good  results  for  all  the  flow  situation. 
Table  2  shows  the  constants  for  all  the  models;  Model  I  refers  to 
the  full  non-linear  expression.  Model  11(1)  and  Model  11(2)  refer 
to  the  the  expression  with  C2  as  a  function  of  the  second  and 
second  +  third  invariant  respectively. 


6.  Results  and  Discussion 


6.1  Return  to  Isotropy 


r 
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Figures  F-01  and  F-02  shows  the  comparison  between 
Uberoi  (1956)  experiments  and  the  predictions  for  the  Model  I  and 
Model  LRR  (  in  this  case  Model  II  reduces  to  LRR  );  it  seems  that 

а.  best  perfomance  is  obtained  when  the  first  part  of  the 
pressure-strain  is  expressed  in  a  non-linear  form.  However  not 
much  improvement  is  gained  for  the  plane  flow  plotted  in  F-03. 

б. 2  Uniform  Plane-strain 

Figures  F-04  shows  the  computed  and  measured  behaviour 
of  the  flow  studied  by  Townsend.  Not  a  great  difference  is 
obseved  in  the  predictions  of  the  three  models  although  Model 
11(1)  shows  better  agreement  with  the  experiments.  Turning  to 
Tucker  and  Reynolds  experiment  figure  F-05  shows  the  perfomance 
of  these  models;  it  appears  that  all  models  fail  to  predict  the 
rise  of  the  v2  component  at  the  end  of  the  contraction  and  in 
particular  the  non-linear  model  whose  rate  of  transfer  of  energy 
among  components  seems  to  be  small. 

6.3  Axisymmetrlc  Contraction 

Figures  F-06,  F-07,  F-08  show  the  computed  and  measured 
results  for  three  cases  taken  from  Tan-atichat  experiments  for 
which  the  strain  rate  is  no  constant.  In  all  cases  studied  the 
models  were  unable  to  reproduce  the  sudden  increase  of  the 
streamwise  component  arising  from  the  acceleration.  For  the  LRR 
model  this  is  consistent  with  the  behaviour  described  by  Launder 
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Reece  and  Rodi  (1975)'  When  observing  Uberoi's  (1956)  16:1 

i* 

axiaynmetrlc  experiment.  For  the  non-linear  model  the  effect  is 
quite  the  contrary*  the  transverse  component  is  overpredicted  and 
the  rate  of  transfer  of  energy  is  small;  the  reason  for  this 
behaviour  seems  to  be  the  introduction  of  the  neu  term  S^  with  a 
negative  value  for  the  respective  constant  and  so  it  makes  a 
negative  constribution  to  u2  and  a  positive  to  v2,  therefore  it 
does  not  tend  to  change  the  trend  of  both  components.  The  same 
situation  is  reproduced  for  Uberoi's  experiments  plotted  in 
figures  F-08  and  F-09.  Figure  F-  13  shows  experimental  results 
and  predictions  using  Model  I  and  Model  11(2)  for  Tan-atichat 
experiments,  it  seems  that  the  introduction  of  the  third 
invariant  has  a  favourable  effect. 

6.9  Shear  turbulence 

Turning  to  Champagne  et  al .  experiment  it  can  be  seen 
in  figures  F— 1 1  that  the  models  reproduce  almost  the  same 
behaviour  but  Model  I  (non-linear)  is  inferior  since  underpredict 
the  shear  component.  Besides  the  level  of  turbulence  energy  at 
the  end  of  the  computation  is  underpredicted  for  all  models. 

Finally  figures  F-12  show  the  shear  flow  of  Harris  et 
al.  ,  here  the  three  models  behave  very  different.  The  u2  and  the 
v2  components  are  overpredicted  by  LRR  Model  ,  Model  I  (non¬ 
linear)  gives  better  agreement  but  model  II  (  C2  variable  )  seems 
to  be  the  best.  Predictions  for  v2  are  unsatisfactory  for  all 
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models.  The  shear  stress  level  is  again  overpredicted  by  LRR 

I 

■odel  and  thii  seems  to  be  the  reason  why  the  predictions  for  the 
other  components  are  unsatisfactory.  An  overall  better  agreement 
is  obtained  with  Model  I  and  11(1)  and  in  particular  in  figure 
F-1H  a  satisfactory  Improvement  is  obseved  when  using  model 
11(2). 

7.  Concluding  Remarks 

A  general  expression  for  the  pressure-strain  term  of 
the  Reynolds  stress  equation  have  been  studied  and  applied  to  a 
number  of  turbulent  hommogeneous  flows.  Considering  possible 
errors  arising  from  the  strongly  coupled  non-linear  differential 
equations  and  the  difficulties  in  optimising  simultaneously  the 
constants  Involved  the  overall  predictions  are  satlfactory.  The 
simpler  version  which  contains  the  invariants  seems  to  be  the 
appropriate  way  for  further  improvements,  this  idea  should  be 
extended  to  the  simpler  model  of  Gibson  and  Launder  (1976)  and 
look  into  the  effects  of  varying  the  "constant"  associated  with 
the  anisotropy  of  the  production. 
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Report  on  the  Modelling  and  Measurements  made  on  a  Puny 

Combustor 

As  detailed  in  the  introduction,  part  of  the  effort  of  the  past  year  has 
been  to  investigate  various  flow  fields  associated  with  an  axial  entry 
type  dump  combustor. 

During  the  past  year  a  dump  combustor  has  been  built  and  is  shown  below. 
This  is  suitable  for  hot  and  cold  flows  and  has  optical  access  at  its 
inlet  and  along  its  1  metre  length.  The  inlet  and  dump  sections  have  80  mm 
and  160  mm  diameters  respectively  and  with  the  air  facility  available, 
inlet  velocities  of  100  M/s  are  possible. 


10  cm 


-C 

1 

Swirl 

Stabiliser 


Dump  combustor  geometry  used  in  the  measurement  study.  The  inlet 
section  is  segmented  allowing  replacement  of  the  swirl  stabiliser 
by  a  baffle  at  the  same  position  or  at  the  dump  plane. 
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Calculations  and  Measurements 

The  work  to  be  described  concerns  measurements  and  calculations  for  the 
following  inlet  conditions  a)  no  stabiliser  b)  a  three  element  baffle  at 
the  dump  plane  c)  a  three  element  bafflr  at  8  cm  upstream  of  the  dump  plane 
d)  a  swirl  stabiliser  at  B  cm  upstream  of  the  dump  plane. 

In  all  cases  k  -  e  modelling  of  the  flow  field  has  been  applied  and  in  the 
swirl  case,  different  levels  of  the  algebraic  model  have  been  used.  In  add¬ 
ition  for  the  baffle  case,  combustion  calculations  have  been  performed. 

The  baffle  used  consisted  of  thin  plate  with  rectangular  obstructions  of 
30Z  total  blockage  area.  For  the  wirier  the  design  was  based  on  the  work 
of  Buckley  et  al  at  AFWAL  and  employed  computed  blade  shapes.  Results  for 
an  approximate  free  vortex  swirler  of  3wirl  number  0.4  are  reported  here. 
The  UDA  system  employed  consisted  of  a  backscatter  arrangement  with  freq¬ 
uency  shifting,  and  signal  processing  was  by  means  of  a  burst  correlator. 
The  light  from  a  488  nm  argon  laser  was  focused  onto  a  rotating  diffraction 
grating  using  a  200  ton  lens.  This  ensures  a  circular  cross  section  of  the 
diffracted  beams  thus  maximising  subsequent  beam  crossing  and  minimising 
uncorrelated  noise.  Two  furhter  lenses  collimate  the  i  orders  and  also 
permit  the  choice  of  beau  separation  and  hence  fringe  spacing. 

The  beams  are  focused  and  crossed  by  means  of  an  output  lens  of  1500  mm 
focal  1  ngth  thus  providing  some  considerable  distance  between  the  test 
rig  and  the  instrumentation.  For  single  photon  collection,  as  with  this 
configuration,  data  processing  is  by  autocorrelation  and  summation,  result¬ 
ing  in  accumulated  autocorrelation  functions  containing  all  velocity  var¬ 
iations  occurring  during  the  sampling  interval. 

Each  sampling  operation  is  controlled  by  a  PET  microcomputer  and  data  is 
analysed  by  a  curve  fitting  procedure. 

The  collection  system,  which  was  aligned  at  15°  to  the  output  axis  as  a 
means  of  obtaining  the  desired  the  longitudinal  resolution,  consisted  of 
a  750  mm  focal  length  collection  lens  of  112  mm  aperture  focusing  onto  a 
collection  pinhole  which  controlled  the  lateral  resolution. 

A  Malvern  Instruments  RF313  photmultiplier,  which  was  capable  of  observing 
single  photons,  was  focused  on  the  collection  pinhole  and  fed  its  output 
to  the  50  nsec  correlator. 

In  a  typical  experimental  run  the  test  rig,  which  is  on  a  motor  driven  tra¬ 
verse,  is  scanned  across  the  region  of  interest,  with  the  fringe  velocity 
being  measured  at  each  sample  point  simultaneously  by  feeding  the  counting 
pulses  from  a  grating  monitor  to  on3  of  the  correlator  storage  channels. 


For  the  present  work,  the  fringe  velocity  employed  was  up  to  130  M/s  dep¬ 
ending  on  the  turbulence  level,  and  the  f tinge  spacing  was  88  y.  The  spatial 
resolution,  which  is  determined  largely  by  the  angle  of  collection  and  waist 
diameter  at  the  probe  yolume,  was  around  1  mm  both  longitudinally  and  later¬ 
ally,  based  On  a  100  times  drop  in  collected  signal  strength  at  these  limits. 
The  air  velocities  used  at  the  80  mm  diameter  inlet  were  between  65  and  75  M/S 
being  limited  by  the  supply,  and  seeding  was  by  means  of  0.5  y  source  particles 
of  titanium  dioxide  generated  in  a  10  cm  diameter  and  120  cm  high  fluidised 
bed. 

Background  noise  limited  the  closest  approach  to  the  walls  to  9.24  mm. 

Results 

Isothermal  Flows 

The  first  results  are  shown  in  Figure  15  and  show  a  comparison  of  calculated 
and  measured  velocities  at  three  axial  stations.  For  each  profile  two  cases 
are  displayed  -  with  the  three  element  at  the  dump  plane  and  with  no  baffle 
stabiliser.  For  the  baffle  case,  the  plane  of  measurement  was  on  a  reflection 
symmetry  plane  at  90°  from  the  mid-line  of  one  of  the  baffles.  Also  included 
in  the  Figure  is  a  comparison  of  the  turbulence  levels  at  one  of  the  stations. 
In  all  casec  the  no  baffle  condition  shows  good  agreement  between  measured 
and  calculated  velocity  and  turbulence,  showing  that  in  this  case,  as  would 
be  expected,  the  k  -  e  model  is  satisfactory. 

For  the  case  with  baffles  at  the  dump  plane,  the  agreement  is  moderately 
good  at  50  mm  from  the  dump  plane  but  deteriorates  as  one  goes  down  stream, 
with  the  measured  values  flattening  out  much  more  quickly.  Examination  of 
the  turbulence  values  indicates  that  the  model  has  greatly  underestimated 
these,  and  most  likely  accounts  for  the  much  slower  levelling  out  of  the 
velocity. 

Comparisons  for  the  swirl  stabilised  case  are  shown  in  Figure  16  and  show 
calculations  based  on  the  k  -  e  turbulence  model  and  the  algebraic  model. 

It  can  be  seen  that  both  models  fall  quite  far  short  of  representing  the 
flow  field  accurately.  However  for  the  algebraic  model,  there  is  a  small 
region  of  reverse  flow  on  the  axis  as  found  experimentally.  This  was  a 
difficult  feature  to  reproduce  and  only  resulted  when  a  full  solution  of 
the  stress  transport  equations  was  applied. 

These  results  contrast  strongly  with  the  simpler  swirling  case  associated 
with  tangential  entry  devices,  which  the  algebraic  model  has  represented 
very  well  and  indicates  that  for  flows  with  swirl  introduced  at  the  axis, 
further  development  work  is  required. 


To  investigate  the  models  further ,  measurements  and  calculations  using  the 
k  -  e  model  have  been  carried  out  on  the  same  geometry  with  the  baffles 
situated  at  8  cm  upstream  of  the  dump  plane. 

Experimentally,  this  jresults  in  the  same  turbulence  levels  but  a  greatly 
reduced  recirculation  bubble  at  the  dump  plane.  Figure  17  shows  the  comp¬ 
arison  results  at  the  150  nm  position  .  It  can  be  seen  here  that  the 
agreement  is  much  better  than  with  the  baffle  at  the  dump  plane  indicating 
that  the  strong  recirculation  in  the  first  case  is  reducing  the  closeness 
of  fit  of  the  k  -  e  calculated  fields. 

Combustion  Case 

Figures  18, 19and  20  show  some  typical  results  which  can  be  obtained  for 
fuel  spray  combustion.  The  three  element  baffle  has  been  used  here  and  a 
single  spray  is  injected  from  upstream  at  right  angles  to  the  wall. 

A  distribution  of  drop  sizes  is  employed  ranging  from  5  y  to  100  y,  with 
a  Rosin  Rammler  distribution  of  mean  50  y  and  width  parameter  of  2.2. 
Preheated  air  at  550°K  is  used  as  input  and  a  fuel  air  ratio  of  0.02  is 
used.  Figure  18  shows  the  trajectory  and  evaporation  history  of  ten  drops 
of  different  size.  The  smallest  penetrate  least  into  the  flow  and  ace 
convested  inside  the  dump.  The  croses  show  the  points  at  which  they  evap¬ 
orate.  It  can  also  be  seen  that  some  of  the  larger  droplets  do  not  evap¬ 
orate  within  the  range  considered  and  it  is  clear  that  this  kind  of  study 
is  of  some  value  in  matching  a  particular  atomiser  to  the  flow  field. 

Figure  19  is  derived  from  the  information  of  the  kind  seen  in  Figure 
for  each  computational  cell  and  shows  a  contour  of  constant  evaporated 
fuel  mass  fraction.  This  is  invaluable  information  in  indicating  where 
unburnt  evaporated  fuel  originates  and  is  convected.  It  can  be  seen  from 
the  figure  that  the  pattern  is  following  the  general  flow  field  and  then 
becoming  entrained  in  the  strong  recirculation  zone  behind  the  baffles. 
Figure  20  is  a  similar  contour  plot  for  temperature  at  2000°K  and  again 
shows  the  curvature  associated  with  one  of  the  recirculation  bubbles.  Such 
information  is  useful  in  association  with  evaporation  patterns  in  correct¬ 
ing  for  hot  wall  regions. 

Conclusions 

Some  of  the  main  features  of  the  flow/combustion  model  have  been  described, 
and  the  ease  of  setting  up  a  device  geometry  has  been  emphasised. 

The  dump  combustor  study  has  underlined  some  of  the  limits  of  applicability 
cf  the  k  -  £  model  and,  in  addition  to  the  axially  symmetric  swirling  flow 


these  appear  also  to  include  flows  with  strong  recirculating  zones.  From 
previous  work,  it  appears  that  a  more  fundamental  approach  to  the  shear 
stress  modelling  may  correct  most  of  the  difficulty  in  these  cases  and 
the  incorporation  of  these  new  developments  for  the  strongly  3-D  case  is 
currently  underway. 

For  the  case  of  swirl  introduced  around  the  axial  region  it  was  seen  that 
the  algebraic  stress  model  of  turbulence  was  able  to  some  degree  reprod¬ 
uce  the  reverse  flow  on  axis,  but  that  there  is  still  considerable  room 
for  improvement  in  this  mathematically  more  difficult  case. 
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TURBULENCE  MODELLING 


FIGURE  1.  OPTIONS  AVAILABLE  FOR  TURBULENCE  MODELLING.1. 
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FIGURE  2  Demonstration  of  system  geometry  specification  within  combustion/ 

flow  algorithm.  This  part  of  the  print  out  shows  the  introduction  of 
a  baffle  at  the  exit  of  a  dump  combustor.  W  -  Wall,  I  -  Inlet. 
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FIGURE  15  Radial  variation  of  axial  velocity  with  and  without  baffle 
stabiliser  situated  at  the  dump  plane.  A)  50  mm  from  the 
dump,  B)  150  mm  from  the  dump,  C)  330  mm  from  the  dump. 

D)  Radial  variation  of  rms  turbulence  at  150  mm  from  the 
dump  plane  with  and  without  baffle. 

O  and  A  are  experimental  points, - are  calculated 

with  coarse  grid,  -  are  calculated  with  fine  grid. 
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FIGURE  16  Radial  variation  of  axial  velocity  at  several  stations  from 
the  dump  plane  with  approximate  free  vortex  swirler. 

© are  experimental  points,  -  are  calculated  using 

k  -  e  model  of  turbulence,  ■  ■  ■  are  calculated  using  the 
algebraic  model  of  turbulence. 
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FIGURE  17  Radial  variation  of  axial  velocity  at  150  ran  from  the  dump  plane 
with  and  without  baffle  stabiliser.  In  this  case  the  baffle  is 
situated  at  8  cm  upstream  of  the  dump  plane. 

□  and  A  are  experimental  points;  - are  calculated  with 

a  coarse  grid, - are  calculated  with  a  fine  grid. 


